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In this paper, we study cooperative scattering of low intensity light by a cloud of N two-level sys- 
tems. We include the incident laser field driving these two-level systems and compute the radiation 
pressure force on the center of mass of the cloud. This signature is of particular interest for exper- 
iments with laser cooled atoms. Including the complex coupling between dipoles in a scalar model 
for dilute clouds of two-level systems, we obtain expression for cooperative scattering forces taking 
into account the collective Lamb shift. We also derive the expression of the radiation pressure force 
on a large cloud of two-level systems from an heuristic approach and show that at lowest driving 
intensities this force is identical for a product and an entangled state. 



I. INTRODUCTION 



Cooperative scattering by an assembly of resonant systems has been studied in detail for many years and is based on 
the seminal work by R. Dicke in 1954 [TJ. Related superradiance effects and collective level shifts have been studied in 
the context of atomic physics in the 70s [1H1| ■ In the last decade, this topic has seen renewed interest [5HT1] with novel 
experiments in nuclear physics |15j and in laser cooled clouds of atoms 16-20J, applications in quantum information 
|21j and quantum phase transitions [22\ 123] . As we are mainly concerned with applications on laser cooled atomic 
samples, we focus in this paper on specific parameters and observables which are of interest in such experiments. We 
therefore derive expressions of the radiation pressure force acting on the center of mass of the atomic cloud, as well 
as the scattered electric field. We go beyond past approximations including the complex kernel for the coupling terms 
between N atoms [21 |5j , described by two-level systems in a scalar approach. Neglecting the complete vectorial nature 
of the dipole dipole coupling seems a priori more justified in a dilute sample of atoms, where near field corrections 
are small [T^]. Furthermore, we obtain the force and the radiation field as quantum operators, which may be useful 
for studying fluctuations and diffusion effects in radiation forces and scattered emission. Also, the imaginary part of 
the complex kernel, describing the collective Lamb shift, is evaluated for a gaussian density profile. 

This paper is organized as follows: in section [TTJ we specify the Hamiltonian used and discuss our approximations. 
In section III we introduce the observables relevant for experiments with cold atoms, namely the radiation pressure 
forces on the center of mass of the atomic cloud and the scattered light intensity. The evaluation of these observables 
is done for specific atomic states in section [TV] We derive the result for this cooperative radiation pressure force from 
a more heuristic approach in section |V| In section |VI| we discuss the relevance of the Timed Dicke State compared to 



a product state for this cooperative pressure force in the low intensity limit before concluding in section VII 



II. HAMILTONIAN AND OPERATOR EQUATIONS 



Our system consists of a gas of N two-level atoms (with random positions Yj, lower and upper states \gj) and \ej) 
with j = 1,...,N, transition frequency uj a with linewidth T = cP : lo 3 / '2i\ he c 3 , where d is the electric dipole matrix 
element), driven by a uniform resonant radiation beam with wave vector ko = k$e z , frequency luq — u a + Ao and 
electric field E (see fig. [I]). 

The atom-field interaction Hamiltonian in the rotating-wave approximation (RWA) is 

= i?o + #i (1) 
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FIG. 1: (color online) Experimental configuration considered : a cloud of two-level atoms is driven by an incident laser detuned 
by Aq from the atomic resonance ui a , with wavevector ko. 
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Here 17 = c?i?o is the pump Rabi frequency, Sk is the photon annihilation operator with wavenumber k and frequency 
ujk = ck, gk = dyjuik/ (2HeoVphj, V p h the photon volume, aj — exp(iAoi)|<?j)(ej| and asj — \ej){ej\ — \gj)(gj\- Instead 
of solving the Schrodinger equation introducing some ansatz for the system state |18j . we write the motion 

equations of the atomic and field operators, 
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We consider the atoms initially in their ground state and we assume weak excitation (f2o <C T), so that we approximate 
63.3(4) ~ I j , where Ij is the identity operator for the jth atom. This approximation amounts to neglect saturation and 
multi-excitation, i.e. all the processes generating more than one photon at the same time (linear regime). Integrating 
Eq. ([5]) and substituting it into Eq. ([3| , neglecting (0) (since the initial field state is vacuum) we obtain 
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~dt 
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The last term in Eq.([6]) describes the effect of the spontaneously emitted photons on the atoms, and it is well known 
in the quantum electrodynamic literature [241 125] . In the Markov approximation (i.e. when the photon transit 
time through the atomic sample is much shorter than the excitation decay time |26|). we assume under the integral 
a m (t — t') ~ a m (t). The time integral then yields a real part (with a term 6(k — ko) ) and an imaginary part 
(corresponding to the principal part of the integral). Taking into account these two terms is at the origin of the 
exponential kernel whereas the real part alone would lead to a sin kernel in Eq.([9| below. We then transform the sum 
over the modes k into an integral, ^ k — ► (Vph/Sn 3 ) J cZk. The real and imaginary parts of the double integral over t 
and k yield the cooperative decay and frequency shift (collective Lamb shift), respectively. The proper expression of 
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the cooperative frequency shift has been obtained adding to the Hamiltonian ^ the not-RWA contributions associated 
to virtual photons exchanged between different atoms. It results the following relation [12] : 
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where T = V ph gl Q kl/ '(ttc). Using Eq.([7b in Eq.([6| we obtain [T2] . 
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and rj m = I i"j — r m | . Eqs.Q describe the time evolution of the atomic operators for N weakly excited atoms scattering 
radiation. The real part of jj m describes the spontaneous emission decay and the imaginary part of 7j TO describes 
the energy shift due to resonant dipole-dipole interactions. A slightly different approach can be used to derive this 
result as shown in appendix [A] Note that even though this result will yield a density dependent collective shift of the 
resonance, we use a scalar model for the held, neglecting thus any polarization and near field dependence [101 113] . 
Detailed calculations for small and large samples of various geometries however show that near field and far field 
contributions as well as resonant and antiresonant terms need to be taken properly into account for quantitative 
predictions [5J [131 [H] , and the present model thus needs to be considered with care illustrating only a part of the 
dipole-dipole coupling for real systems. 
Eq.([8| can also cast in the form 

<la > 1 : {a J ,H>+H eff }, (10) 
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and the commutation rules in the linear regime are [&j , crjj = S 



III. OBSERVABLES 

Among the different observables of the system, scattered light and radiation pressure force contain important 
signatures of cooperative scattering. Concerning scattered radiation, the positive-frequency part of the electric field 
is defined as 

E(r, t)=iJ2 ^fik(i)e-^ t+!k ' r (12) 



where £k — ^JhuJk/2eoV p h is the single-photon electric field. By integrating Eq.([5j) and inserting it in Eq.(12) we 
obtain 

E(r, t) = J2 £k9k Y\ e ik < r - rm) - iuot / dt' er 1 ^-^ 1 ' a m (t - t') (13) 
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Using Eq.([7|), the Markov approximation leads to 
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which has a transparent interpretation as the sum of wavelets scattered by N dipoles of position Vj and detected at 
distance r and time t. In the far field limit, |r — r } | r — (r ■ Yj)/r and 
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where k s = fco(r/r). 

The radiation pressure force acting on the jth-atom has been calculated from Eq.|l| as Fj = ~V Tj H = F a j + F e j 
where ITS] 
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where F aj - and F ej - result from the recoil received upon absorption of a photon from the pump and from the emission 
of a photon into any direction k, respectively. Eliminating the field using Eq.([5]), Eq.(17) becomes 
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Assuming the Markov approximation, a m (t — t') ~ cr m (t), then Eq.(18) becomes 
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where r jm = rj — r m . The force (19 1 acting on the jth atom has a single-atom contribution (term m = j in the 

sum) accounting for its own photon emission recoil, and a contribution F*,^ (terms m ^ j) accounting for coupling 
between the jth. atom and all the other atoms. Note that this dipole-dipole interaction can occur via a coupling to 
common vacuum modes of radiation. The interference terms in the total scattered field can leave a fingerprint on the 
forces acting on the atoms inside the cloud. The first contribution yields 
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where the sum is over all the randomly oriented modes k — k k and we have omitted the self-energy shift (Lamb shift) 
coming from the principal part term of the time integral in Eq( 19 ). Noting that for m ^ j we have zk exp(zk • r jm ) = 
V rj exp(ik • fjm), the second contribution to Eq.(19) can be written as 
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Using Eq.Q in Eq.Q, Eq.([l9]) becomes 
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is the effective interaction energy between jth and mth atoms. Since V r [exp(zfcor)/r] = r(ikor — 1 ) exp(i/cor)/r 3 , 
Eq.(22 1 becomes 
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where iij m — Yj m / r jm- The emission force acting on the jth atom has two contributions: a) a self-force due to its 
own photon emission; b) a force due to the dipole-dipole interactions with all the other atoms. This second force has 
a term decreasing as l/rj m and one decreasing as l/r jm 2 . 



IV. ATOMIC STATE 



The linear approximation assumed in the equations of the atomic operators <Xj, Eq.(|8|), suggests that we may 
restrict the Hilbert space of the N atoms to the subspace spanned by the ground state \g) = \g%, . . . , g^) and the 
single-excited-atom states \j) — \gi, . . . , ej, . . . , ^at) with j = 1, . . . , N. Hence, we set 
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where we will approximate a ~ 1 after having evaluated the different expectation values, e.g. (aj) sa /3j and 



(<rja m ) « (3*p m . So, Eq.^ yields 
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with initial conditions j3j(0) — 0. The self-interaction term, Tjjj = T — iAQ^s yields the single-atom spontaneous 
decay F and the single-atom Lamb shift Ail^s, which can be reabsorbed in the definition of the atomic frequency ui a , 
and will be neglected in the present analysis. 



Considering the force applied to the center-of mass of the atomic ensemble, F = (1/iV) ■ Fj, from Eqs.(16) and 



( 24 1 the components along the z axis are 
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where ji(z) = sin(z)/z 2 — cos(z)/z is the first order spherical Bessel function and Zj 
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(28) 
Note also 



that the self-force (|20[) has zero average since J2 k k = (although in general its fluctuations are different from zero). 

Also, from Eq.(|l5| it is possible to obtain the average intensity of the scattered radiation as a function of the atomic 
wave function, 
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The state (25) may be conveniently expressed in the timed Dicke (TD) basis, introduced originally by Dicke pQ and 
successively considered by R. Friedberg and coworkers [3) for their study on cooperative Lamb shift and by M.O. 
Scully and coworkers [9] to describe cooperative decay of TV atoms prepar ed i n a symmetric phased state. The 
completely symmetric TD state is |+)k = J^.- exp(iko • r j ) | j ) and Eq.(25) can be written as 
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A numerical analysis of Eq.(26) shows that, for a constant driving field Qq and for atomic cloud sizes much larger 
than the optical wavelength, the occupation probability of the states |s)k is only a small fraction of the atomic state 



and it is in general negligible, so that Eq.(26) becomes 
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Lamb frequency shift [31 [TO] . At steady state we find 



is the factor form and the integral over k in Eq.(33) is evaluated as a principal part. The term A^v is the collective 
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The cooperative radiation force can be obtained from the standard single-atom radiation pressure force Fi = 
hk rQQ/(4:A^ + T 2 ) substituting the natural linewidth by the collective linewidth, r^v = TNsn, and multiplying it 
by 1 — /n/sn , where /n / sn is the probability to observe the photon emitted in the forward direction. Isolating the 
term j = m, 



\S N (k)f = 1 + e^"-^-^ « 1 + ^(k) 
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where the factor form Soo(\s) is evaluated for a continuous approximation with density distribution n(r), 

Soo(k) = i y dm(r)e J ( k »- k ) r . (39) 
Then, ps (l/N) + Soo and ~ foo where, 
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FIG. 2: (color online) Distributions for values of sjv, /jv and Ajv for N = 50 atoms, plotted for 10000 configurations for a size 
corresponding to a — 10 (blue curves) and a — 2 (green curves). 
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FIG. 3: (color online) Collective Lamb shift vs atom number for a = 1.6 (green triangles) a — 2 (blue circles) and a = 3 (red 
squares) . 



The factor form SoaQ*) and the integrated factors Soo and foo have been calculated in ref. |18j for a Gaussian 
density distribution with ellipsoidal profile, n(r) exp[—(x 2 +y 2 )/2a 2 — z 2 /2a 2 ], yielding S 00 (ko, 9) = exp{— a 2 [sin 2 8 + 
rj 2 (cos9 — l) 2 ]/2}, where a = k a r and T] = cr z /a r is the aspect ratio. For spherical and large clouds (77 = 1 and 



B) 



l/(4a 2 ), f a 



l/(8cr 4 ) and the collective Lamb shift is Ajv ~ A-oo where (see [Hj and Appendix 



FN 
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which is a redshift, proportional to the number of atoms in a cubic wavelength [3], i.e. atomic density and not optical 
thickness bo — 3N/a 2 . These values for s,^, f x and can be compared to numerical evaluation of the sjv, /jv and 
Ajv for a finite number of atoms and a specific configuration. In Fig. [2] we show the distribution of these values for 
different sample size. 

In our numerical simulations shown in Fig. [3] we observe strong configuration dependent fluctuations for the value 
of the collective Lamb shift. A precise comparison with our analytical expression, valid for large clouds, is thus 
cumbersome and did not allow us to validate precise predictions of the numerical factor in Eq. ( 42 ) . 



Normalizing the radiation pressure force with respect to the single atom value, we obtain for large atomic samples, 
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FIG. 4: (color online) Emission diagram computed according to Eq. (29 1 for the Timed Dicke state |+}k with N = 40 atoms 
: a — 0.4 (blue), a = 1 (red) , a = 8 (green). 



Finally, from Eq.(29) we obtain the scattered intensity 
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This expression of the scattered intensity illustrates the role of the shape of the atomic cloud for the modified emission 
diagram. The emission diagram of the TD state is shown in Fig. [4] It illustrates the strong forward emission by the 
cloud when its size exceeds a few optical wavelengths, reminiscent of Mie scattering, or more precisely of Raylcigh- 
Debye-Gans [28] . As we will discuss in the following section, a modified emission diagram yields a modified radiation 
pressure force, as the recoil of the scattered photon (partially) compensate the recoil effect at absorption. 



HEURISTIC APPROACH 



The result ( 36 ) can be interpreted heuristically considering the momentum balance in a given time interval St |29j . 
During St, N two-level atoms with positions r 3 - (j = 1, . . . , N) do SN florescence cycles, each time absorbing a photon 
with momentum Kk from the laser and emitting a photon with momentum filq [i = 1, . . . , SN) in a random direction 
fij, with probability Pjj = P(fli,rj). The momentum variation for the jth atom after SN cycles is 



SN 



S Pj = {hk )SN - Y,{hki)P, 



h3- 



(45) 



For a single isolated atom the emission is isotropic and Pjj = 1, but for N atoms the emission can be not isotropic 
depending on the atomic distribution. Also, the excitation could be not uniform if the phase front of the driving beam 
is getting distorted by the refractive index changes in the atomic cloud. Assuming for simplicity that the excitation 
is uniform over the entire atomic ensemble and neglecting phase distortion effects [201 [HI] , SN will be the same for 
all the atoms and |kj| = kg. Considering the momentum variation along the direction of the incident photon (z axis), 
after averaging over the atoms 



SN 
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where Pj = (1/N) . Pjj = P(cos#i) is the emission probability along the angle 0t. Considering cos^i and SN as 



independent random variables, the statistical average of Eq.(46) is 
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where we assumed J2i cos ®i ~ $N ■ cos 9. Hence, the pressure force is 
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Comparing with Eq. ( 36 1 we found the following correspondence 
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where — TNsn- So, the scattering rate (8N/St) is equal to the excitation probability, p ee — f2g/[4(A — Aat) 2 - 
times the collective decay rate, T^. The radiation pressure force (36 1 is equal to the momentum photon, fi,fc n , multiplied 
by the scattering rate and by a geometrical factor 1 — cos 9 taking into account the directionality of the scattered light. 
Cooperativity modifies both the scattering rate, enhancing the decay rate and shifting the resonance frequency, and the 
scattering direction. Small samples tend to radiate isotropically whereas extended samples radiate superradiantly in 
forward direction [111 I18j . These cooperative effects can be revealed measuring radiation pressure force by monitoring 
center-of-mass motion of large atomic clouds released by magneto-optical traps [TH [17], and then identifying fast 
decay, shifts and modified emission diagrams described by Eqs.(36) and (44 1. 



VI. PRODUCT STATE 



It has been noted that the same results obtained for a symmetric TD state could be obtained assuming a product 
state for N atoms 0[T3] (named also 'Bloch state' by some authors [7]): 



\^(t)) c = H{a c (t)\g J )+(5 c (t)e 



iko-r,- 



(50) 



where a c (t) and ft c (t) are the same for every atom, with |a c (i)| 2 + |/3 c (i)| 2 = 1. The ansatz of Eq. (50 1 assumes each 
jth atom driven into the excited state with equal probability |/3 C (<)| 2 and phase (j)j = ko • Vj — Aoi. As it happens 



for the symmetric TD state (30), the driving field imposes a coherence in the photons emitted spontaneously by each 
atom, so that superradiance arises because the state is symmetric under exchange of particles [32] . However, it is 
expected that the quantum statistic of the symmetric TD state will be quite different from that of the 'quasi-classical' 
product state. Notice that for |/3 C | -C 1 the product state (50) can be written in the following form [71 IT4] 
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where \j, m) = \g\, . . . ,ej, . . . , e m , . . . ,gjsr). Hence, the product state can be expanded in the symmetric TD states 
with 1 to TV excited atoms. Only in the limits a c « 1 and |/3 C | <C 1 the product state reduces to the symmetric single- 
excited atom state \tp) « \g) + /3 c -\A/V|+)ko if only the first two terms of Eq. (51) are retained. The expectation values 
for the state (50) are (oj) = a*/3 c and (a^&j) = \a c \ 2 \f3 c \ 2 , so for a c s=s 1 they coincide with those obtained from the 



symmetric TD state. Differences between the product and the symmetric TD states should appear when higher-order 
expectation values are observed, as for instance (&ja m ), which is zero for the TD state and a* 2 /3jf3 m sa f3j(3 m for the 
product state. Notice that operator ordering produces different results in high-order expectation values if scattered 
photons or atomic forces are measured. These features and non classical effects studies in cooperative scattering by 
cold atoms will be the object of a future investigation. 



VII. CONCLUSION 



In this paper, we have included a more precise kernel to evaluate the cooperative radiation pressure force on a 
cloud of two- level systems. The collective Lamb shift leads to a shift Ajy of the resonance, which is proportional to 
the spatial density. As we have used a scalar model in this paper, near field and polarization effects are neglected. 
One thus needs to consider this shift with some scepticism as the numerical factor for this shift in a real system 
will be strongly modified by the vectorial nature of the light For dilute clouds, we recover previous results |16j . 
where these density effects are negligible. We also presented a simple model to estimate the radiation pressure force 
from the modified emission diagram and assuming coupling to the single photon superradiant (Timed Dicke) state 
[5] . This approach can be useful to estimate not only average forces but also fluctuations and dissipation. Finally, we 
noted that in the low intensity limit, the average result we derived for the cooperative radiation pressure force can be 
obtained either by assuming a driven Timed Dicke state or a product state [5) ITj IT3]. with no entanglement required. 
Looking for non classical features in cooperative scattering of light by a cloud of two-level system thus requires studies 
of higher orders either by using higher intensities or looking at correlations or fluctuations of the force. 
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Appendix A: Evaluation of the integral kernel in Eq.([6]) 

Let's consider the last term in Eq.Q and pass to the continuous frequency approximation: 

I(r iro ) = 5>^ ik - r ^ J^dt'a^t-qe-^-^' ^ 
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We exchange the integration order and introduce spherical coordinates, dk = dkk 2 d<f> d9 sin 9. After integration of 
the angular part, we obtain 
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where r jm = |rj m |. We approximate the k integral as 
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where we made the following approximations: a) we assumed the spectrum centered around k m kg, so that kg 2 ~ 
kog 2 g ; b) we extended the lower integration value from to — oo, since the relevant values of k are around kg. Using 
the expression above, we write 



I( T jm) 



dt'a m (t - t')e ickot ' {6(f - r jm /c) - 6(t> + r jrn /c)} = 



korj 



T e ltc ° 
2 ik^rj 



-a m (t~r jm /c). (A4) 



where T — Vphk^g 2 . / (7rc). We observe that this approach does not require to assume the Markov approximation 
before solving the time integral, as in the standard approach |12j . On the contrary, this approach allows to obtain the 
retarded (or not local) kernel, which, when the 'rapid transit approximation' is assumed, i.e. a m (t — rj m /c) s» a m (t), 
reduces to the exponential kernel of Eq. (8) . 



Appendix B: Collective Lamb shift for a Gaussian distribution 



Let consider Eq.(33) for a continuous distribution: 

rjV f°° dun 3 '** 



8tt 2 J k-IJo 



d6 sin6\Soc(K,0,' 



(Bl) 



A spherical Gaussian distrib utio n, n(r) exp(— r 2 /2aj t ), yields S^k, 9, 0) = cxp[— <j 2 (k 2 + 1 — 2kcos0)/2], where 



a = ko<Tn. Inserting it in eq.(Bl I we obtain 



< ' tJ(;J j ) r d9 sin 9e 2 
o 



47T Jq K — 1 

TN r°° dun 2 

^o 1 J 

TN f°° , ( I 
- — / dn « + l 



o 



K — 1 



TN [°° / 1 
r P / dx 2 + x + - 



8na 2 



e -a\^if _ e -<r 2 ( K +i) 2 



-a 2 (2+x) 2 



(B2) 
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For a 3> 1 it is approximated by 



in agreement with the result of Friedberg and Manassah [Tl] . 



[1] R. H. Dicke, Phys. Rev. 93, 99 (1954). 
[2] R. Lehmberg, Phys. Rev. A 2, 883 (1968). 

[3] R. Friedberg, S.R. Hartmann, J.T. Manassah, Phys. Rep. 7, 101 (1973). 

[4] M. Gross, S. Haroche, Phys. Rep. 93, 301 (1982). 

[5] J. H. Eberly, J. Phys. B: At. Mol. Opt. Phys. 39, S599 (2006). 

[6] M. O. Scully, E. S. Fry, C. H. R. Ooi, and K. Wodkiewicz, Phys. Rev. Lett. 96, 010501 (2006). 
[7] R. Friedberg, and J.T. Manassah, Laser Phys. Lett. 4,900 (2007). 

[8] A. A. Svidzinsky, J.-T. Chang, and M. O. Scully, Phys. Rev. Lett. 100, 160504 (2008). 

[9] M. O. Scully and A. A. Svidzinsky, Phys. Rev. Lett. 373, 1283 (2009). 
[10] M. O. Scully, Phys. Rev. Lett. 102, 143601 (2009). 
[11] S. Prasad, R.J. Glauber, Phys. Rev. A 82 063805 (2010). 

[12] A. A. Svidzinsky, J.-T. Chang, and M. O. Scully, Phys. Rev. A 81 (2010) 053821. 
[13] R. Friedberg, and J.T. Manassah, Phys. Rev. A 81, 063822 (2010). 
[14] R. Friedberg, and J.T. Manassah, Phys. Lett. A 374 1648 (2010). 
[15] R. Roehlsberger et al., Science 328, 1239 (2010). 

[16] T. Bienaime, S. Bux, E. Lucioni, Ph. W. Courteille, N. Piovella, and R. Kaiser, Phys. Rev. Lett. 104 (2010) 183602. 
[17] H. Bender, C. Stehle, S. Slama, R. Kaiser, N. Piovella, C. Zimmermann, and Ph. W. Courteille, Phys. Rev. A 82 (2010) 
011404. 

[18] Ph.W. Courteille, S. Bux, E. Lucioni, K. Lauber, T. Bienaime, R. Kaiser, N. Piovella, Eur. J. Phys. D 58, 69 (2010). 
[19] R. Kaiser, J.Mod. Opt. 56, 2082 (2009). 
[20] S. Bux et al, J.Mod. Opt. 57, 1841 (2010). 

[21] A. Greeentree, C. Tahan, J. Cole, L. Hollenberg, Nat. Phys. 2, 856 (2006). 
[22] A. Osterloh, L. Amico, G. Falci, R. Fazio, Nat. 416, 608 (2002). 
[23] E. Akkermanns, A. Gero, R. Kaiser, Phys. Rev. Lett. 101, 103602 (2008). 
[24] M.O. Scully and S. Zubairy, Quantum Optics, Cambridge Univ. Press, 1997. 

[25] G.S. Agarwal, Quantum Statistical Theories of Spontaneous Emission and their Realation to other Approaches, Spriger 
tract in Modern Physics, ed. G. Hohler, Springer- Verlag, Berlin 1974. 

[26] Light may propagate in dense atomic samples with a group velocity smaller than c 27 . In these case the Markov approx- 
imation should be satisfied by a more stringent condition. 

[27] G. Labeyrie et al., Phys. Rev. Lett. 91, 223904 (2003); 

[28] H. C. van de Hulst, "Light Scattering by Small Particles", Dover Publications Inc., New York (1981). 

[29] J. Dalibard, PhD, Universite Pierre et Marie Curie - Paris VI (1986). 

[30] J. Gordon, Phys.Rev. A 8, 14 (1973). 

[31] G. Campbell et al., Phys.Rev. Lett. 94, 170403 (2005). 

[32] M. Sargent III, M.O. Scully, W.E. Lamb, Laser Physics, Addison-Wesley Publ. 1974, p.400. 



